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ABSTRACT 


The  basic  definitions  for  viscoelastic  stress  intensity  factors  are  set 
forth  with  associated  computations  for  determining  them  in  laboratory  specimens 
Laboratory  tests  have  been  conducted  to  characterize  propellants  fracturewise, 
and  the  most  suitable  test  found  was  the  long  strip  subjected  to  a  constant 
force  loading.  The  influence  of  dewetting  on  the  fracture  characterization 
appears  to  be  significant,  but  the  matter  has  not  yet  been  fully  resolved. 

An  approximate  method  for  performing  viscoelastic  fracture  analyses  has 
been  presented  which  incorporates  the  fore -mentioned  fracture  characteriza- 
tion. It  requires  minimum  computer  usage  and  should  yield  good  results. 
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VISCOELASTIC  ANALYSIS  METHOD  FOR 
ROCKET  MOTORS  CONTAINING  CRACKS 
INTRODUCTION 

Solid  propellant  has  "been  found  to  "be  susceptible  to  cracking  during 
fabrication  or  service  just  like  other  engineering  materials.  The  cracks 
may  or  may  not  be  detrimental  to  the  rocket  motor  depending  upon  location, 
size,  and  propagation  rate;  however  a  potential  danger  always  exists  when 
a  crack  is  present;  for  instance,  the  crack  introduces  additional  burning 
which  may  cause  significant  overpressure  and  a  subsequent  malfunction.   It 
is  also  possible  that  a  crack  may  create  a  pathway  for  burning  to  the  case 
prematurely,  producing  a  burn-through.   On  the  other  hand,  it  may  lead  to 
other  areas  of  the  grain  where  malfunctioning  occurs  due  to  alterations  in 
the  ballistics  or  due  to  loss  of  structural  integrity  resulting  from  the 
high  stresses  and  strains  associated  with  the  crack  front  exceeding  the 
design  limits.  The  consequences  of  such  developments  must  be  determined 
for  each  individual  situation;  therefore  a  theoretical  basis  for  predicting 
the  initiation  and  propagation  behavior  of  cracks  in  propellant  is  needed. 

The  elastic  theory  of  fracture  mechanics  has  been  successfully  applied 
to  rocket  motor  geometries  and  loadings  by  Francis,  Lindsey  and  Parmerter  (l) 
Similar  work  was  done  on  viscoelastic  materials  by  Francis  et  al  (2)  using 
an  unfilled  polyurethane  for  the  grain  material.  Two-dimensional  grain 
cross -sectional  slices  containing  cracks  were  analyzed  and  compared  with 
experiment.  Predictions  of  viscoelastic  crack  behavior  were  made  by  using 
a  sequence  of  elastic  crack  propagation  problems.  The  method  of  fracture 
characterization  followed  that  used  by  Knauss(3)  on  the  same  polyurethane 
material.     Progress  was  made  and  results  of  the  study  were  used  to 
provide  guidelines  for  some  of  the  concepts  developed  in  this  work. 


In  this  report  the  fundamental  expressions  for  stress  and  displacement 
near  a  crack  tip  are  developed  for  linear,  unfilled,  viscoelastic  materials. 
These  expressions  are  used  to  reduce  crack  propagation  data  obtained  on  strips 
of  propellant  that  are  considered  to  be  infinite  in  extent.  Both  constant 
stress  and  constant  displacement  boundary  conditions  are  employed  in  the  tests 
for  purposes  of  comparison.  Experimental  fracture  characterization  in  terms 
of  stress  intensity  factors  versus  crack  velocity  in  propellant  is  the  end 
result.   Finally  a  method  is  presented  for  the  viscoelastic  analysis  of  rocket 
motors  containing  initial  cracks. 


DEFINITIONS  AND  ASSUMPTIONS 
Stress  Intensity  Factor 

Elastic  solutions  for  the  stress  field  in  the  vicinity  of  a  crack  by 
Williams   ,  Westergaarcr  '  and  Irwin  '   have  been  written  in  terms  of 
undetermined  coefficients  that  became  explicit  when  the  particular  body 
geometry  and  loading  are  specified.   In  terms  of  the  mixed  coordinates 
(see  Fig.  l) 


Figure  1  Stationary  crack  coordinate  system. 


normally  used  to  express  the  stress  field, 


CT,,(r,e)  =  — z=Z  fi^   +  — —   Si^9) 
1J       J  2r   1J      J~2T     1J 


(1) 


where  K,  and  Kp  are  the  opening  mode  and  shearing  mode  stress  intensity 
factors  (S.I.F.)  respectively. 

To  more  precisely  define  these  quantities,  use  can  be  made  of  the 
opening  and  shearing  stresses  in  expanded  form. 
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(2b) 


These  two  stress  components  possess  a  certain  assymetry,  and  the  S.I.F.  can 
he  conveniently  defined  in  terms  of  them. 


K  =  lim  v^r  a  (r,0) 
r  -  0      y 


K_  =  lim  J&   t   (r,0) 


r  -   0 


xy 


(3a) 


(3b) 


where  the  limit  is  used  to  indicate  that  this  relationship  only  applies  in 
a  small  neighborhood  of  the  crack-tip. 

For  6  =  0,  r  =  xf  and  y  =  0.  Making  this  substitution,  the  definition 
can  be  written  also  in  terms  of  cartesian  coordinates. 


K..  =     lim  J2ri  a   (x',0) 


K_  =     lim  J2tf  t   (x',0) 


x  -  0 


xy 


(Ua) 


In  viscoelastic  fracture  no  standard  definitions  of  stress  intensity 
factor,  fracture  toughness,  coordinate  orientation,  etc.  have  emerged; 
however  it  appears  straight  forward  to  extend  the  above  definitions  in  the 
obvious  way.  Since  the  crack  is  normally  propagating  under  load  in  visco- 
elastic materials,  it  is  convenient  to  shift  the  base  of  reference  to  the  system 
shown  in  Fig.  2.      y 


c(t) 


Figure  2  Moving  crack  coordinate  system. 
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Translating  the  previous  expressions  for  S.I.F.  "by  an  amount  c(t)  and 
allowing  for  the  possibility  of  the  stresses  varying  with  time, 

Kv(t)=   lim    £[2(x  -  c(t))]2  a  (x,0,t)}  (5a) 

x  -  c(t)  y 

K  V(t)  =   lim    {[2(x  -  c(t))]2  T   (x,0,t)}  (5b) 

x  -  c(t)  ^ 

Apart  from  the  square -root- of -two  factor,  this  definition  is  the  same 
as  used  by  Grahanr   .  Knauss^  '  also  used  this  form,  but  Meuller^  *   and 
Knauss  and  Meuller    used  a  nondimensional  form  of  the  above  for  an 
infinite  strip  of  width  2b. 

I  =   lim    j[x  *{t>   2  -£ \  (6) 

x  -*  c(t)  LL    b    J     CTo    J 

Since  all  of  these  forms  are  closely  related,  and  since  the  definitions 
of  equation  (5)  possess  the  appropriate  limit  check  for  elastic  materials, 
they  will  be  used  throughout  this  report. 
Fracture  Postulate 

Fracture  studies  in  recent  years  have  led  to  the  conclusion  that  flaws 
in  materials  will  grow  when  under  load.  This  continues  until  the  flaw 
becomes  of  critical  size,  at  which  time  the  specimen  undergoes  unstable 
crack  propagation  and  fracture.  This  has  been  found  to  be  true  for  metals, 
especially  in  fatigue,  and  for  nonmetals  alike. 

Knauss    has  published  theoretical  developments  of  these  ideas  for 
viscoelastic  materials  and  has  shown  that  the  important  parameter  of 
interest  for  rate  sensitive  materials  is  crack  velocity.   In  the  elastic 
case,  the  threshold  of  rapid  crack  propagation  is  often  the  focus  of 
interest;  however  in  viscoelastic  materials  the  crack  is  found  to  propagate 


for  almost  any  load  level.   The  velocities  are  infinites imally  small  for 
low  loads  and  increase  as  the  body  is  more  highly  stressed.  The  question 
to  be  answered  then  becomes,  not  at  what  load  level  does  the  crack  propa- 
gate, but  at  what  velocity  does  the  crack  propagate  for  any  given  load. 
As  defined,  the  viscoelastic  stress  intensity  factor  contains  information 
relative  to  the  specimen  geometry,  material  and  boundary  conditions; 
therefore  we  desire  the  relationship  between  stress  intensity  factor  and 
crack  velocity. 

This  brings  us  to  the  fracture  postulate: 

If  t^^Body  1  "  ^^Body  2 

then  f^W  1  =  [c(t)]Body  2 

In  words  this  states  that  for  every  stress  intensity  factor  there  is  a 
corresponding  unique  crack  velocity. 

Knauss  J   has  performed  experiments  on  Solithane  113?  a  polyurethane 
elastomer,  which  have  substantiated  this  to  a  degree  (Francis  et  al 
have  duplicated  these  tests  and  done  additional  work  on  other  types  of 
loads.)  For  all  temperature  ranges  tested,  the  material  was  elastic  in  the 
sense  that  it  did  not  relax;  however  the  crack  behavior  was  viscoelastic 
and  exhibited  a  velocity  that  varied  with  loading.   Knauss  postulates  that 
in  this  type  of  material,  when  the  strain  at  the  crack  tip  is  the  same  in 
two  different  bodies,  the  same  crack  velocity  will  be  produced.  Of  course 
in  a  body  that  is  behaving  elastically  so  far  as  stress  and  strain  are 
concerned,  the  same  strain  level  is  equivalent  to  the  same  stress  intensity 
factor.  Thus,  in  the  postulate  of  fracture,  we  are  following  the  previous 
ideas  of  Knauss.  Furthermore,  it  states  that  if  the  stress  intensity  factor 
history  can  be  calculated,  the  crack  velocity  equation  can  be  determined 


from  the  experimental  characterization  of  the  material.   Once  this  is 
known,  the  crack  history  for  a  rocket  motor  can  "be  obtained.  This  will 
be  demonstrated  later. 

EXACT  SOLUTIONS  FOR  S.I.F. 

The  classical  correspondence  principle  normally  applied  in  the  solu- 
tion of  boundary  value  problems  in  the  linear,  quasi-static  theory  of 
viscoelasticity ,  does  not  apply  for  crack  propagation  problems.   In  this 
instance  new  boundary  is  being  created  as  the  crack  propagates ,  and  the 
conditions  on  the  boundary  are  changing  in  type;  i.e.  at  the  crack  tips 
they  are  changing  from  displacement  to  stress  boundary  conditions  and  vice 
versa.   In  this  kind  of  problem  there  will  be  points  of  the  boundary  where 
only  partial  histories  of  the  boundary  conditions  are  known  a  priori; 
consequently  the  integral  Laplace  Transform  cannot  be  computed,  and  the 
correspondence  principle  is  not  applicable. 

Graham     has  extended  the  correspondence  principle  to  encompass 
crack  problems .  Two  examples  that  are  of  particular  interest  in  fracture 
characterization  of  propellant  are  the  long  strip  with  either  displacement 
or  stress  boundary  conditions. 
Long  Strip  with  Stress  Boundary  Conditions 

The  sheet  and  loading  are  symmetric  about  the  x  axis,  as  shown  in 
Fig.  3,  where  the  lower  half  of  the  stip  is  removed  and  replaced  by  the 
stresses  which  act  along  the  sheet's  horizontal  centerline. 
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rr~n 
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l  _t 
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A(i) 


Fig   3.      Constant    Siren,     Sheet 


Let  the  entire  boundary  of  the  sheet  be  defined  as  £  with  $0{t) 
representing  the  uncracked  portion  of  the  x-axis  and  p,(t)  representing 
the  remaining  boundary.  Then  £  =  p,(t)  U  Pp(t).  The  boundary  conditions 


are: 

a    (X,t)  =  0  on  p 

s 

(0  on  crack  edge 

an(X,t)  =<0  on  x  =  0,L    ^  on  p^t)  (7) 

a&(X,t)  on  y  =  b 
un(X,t)  =0  on  p2(t) 

the  subscripts  s  and  n  respectively  refer  to  normal  and  shear  components. 
The  boundary  conditions  of  Equation  (7)  meet  the  criteria  of  Graham's 
extended  correspondence  principle.  The  extended  correspondence  principle 
permits  the  construction  of  the  viscoelastic  solution  for  a  moving  crack 
from  an  elastic  solution  for  a  moving  crack. 

CTJ(t)  =  if1  [a.f * (s)]  E  -.  E*  (8) 

u  -*  u* 

where  m  designates  a  moving  crack.   The  moving  crack  solution  may  be 
generated  from  a  similar  stationary  elastic  crack  solution.  For  an 
incompressible  elastic  strip,  the  stresses  in  the  neighborhood  of  the  crack 
tip  are  given  by, 

e    sLiZE_f  .-(e)  (9) 

1J   V2rrr   1J 

where  o  =   boundary  stress  in  uniform  stress  region 

2b  =  strip  width 

The  crack  is  now  imagined  to  propagate  slowly.  For  most  elastomers 
and  all  propellants,  the  velocity  of  propagation  is  quite  low,  which  makes 
the  inert ial  terms  in  the  equations  of  motion  small  and  negligible.  The 
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field  equations  become  the  same  as  those  of  the  static  crack  problem; 
therefore  the  stress  field  is  the  same  as  that  given  in  Equation  (9)? 
except  that  the  coordinate  system  to  which  the  above  quantities  are 
referenced  is  now  moving  since  it  is  attached  to  the  crack  tip.  The 
dynamic,  elastic  solution  is 

K  (t)  ■jfeyE.f  re(t)i  do 

1J  V2wr(t)      1J 

It  is  a  simple  matter  to  reference  these  quantities  to  an  inertial 
coordinate  system.  As  stated  previously  the  loading  and  geometry  is  taken 
such  that  only  the  opening  mode  of  propagation  is  present;  therefore  from 
Fig.  2 

r(t)  cos  9  (t)  +  c(t)  =x  (11) 

Equation  (10)  for  the  dynamic  elastic  problem  becomes , 

me   }  =   a(t)Vb-        (t)]  (12) 

1J      ^2n(x-c(t)   1J 

To  use  the  correspondence  principle,  the  transform  of  Equation   (12)    is 
needed. 


me*  _  r       2&ME        p.  .[6(t)]   1  (13) 


a 


By  Equation  (8),  the  corresponding  viscoelastic  expression  is 


rv*  =  r  a(t)/b"   _  p>  [e(t)]l  (iU) 

1J   V2n[x-c(t)]   1J      J  E  -  E* 


Since  E  and  v  do  not  appear,  the  inversion  is  trivial,  and  the  viscoelastic 
stresses  are  identical  to  the  elastic  ones. 


aV.(t)  = 


g^Wb   F..fe(t)i 


V2rr(x-c(t))   1J 


(15) 


From  this  result  and  by  application  of  the  definition  of  stress  intensity- 
factor  given  in  Equation  (5), 


K^(t)  =  a(t)J  -n 


(16) 


This  viscoelastic  solution  is  accurate  to  the  degree  that  the  elastic 


solution  from  which  it  was  constructed  is  accurate.   In  this  case  it  is  known 

(12) 
to  "be  good  via  Knauss      series  solution. 

Long  Strip  with  Displacement  Boundary  Conditions 

Many  elastic  crack  propagation  studies  have  been  conducted  on  strips 
with  displacement  boundary  conditions,  but  in  order  to  use  it  for  viscoelastic 
fracture,  the  viscoelastic  solution  is  needed.   This  is  not  directly 
obtainable  from  Graham,  but  it  is  possible  to  gem  ate  it  by  following  some  of 
Graham's  methods. 

Consider  the  viscoelastic  boundary  value  problem  of  Fig.  k,   where  A,  B, 
and  C  do  not  involve  material  constants,  and  the  crack  is  propagating. 


uV  =  0 
n 


u 


-  A(x,t)   v  _ 


V 


=  0 


uV  =  B(x,t) 
n 


uV  =  0 
s 


Crack  Surface 


uV  =  0 
n 


s 


V 


u  =  C(x,t) 


Fig.  k     Viscoelastic  problem. 
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Following  Graham,  we  consider  a  dynamic  elastic  problem  with  a  moving  crack 
and  the  same  boundary  conditions.   This  is  shown  in  Fig.  5- 


ue  =0 
n 


u  =  A(x,t) 


u  =  B(x,t) 

n      ' 


u6  =  0 

S 


a?,(x,t) 


u±(x,t) 


ue  =  0 
n 


ae  =  0 

S 


u  =  0 


u  =  C(x,t) 


Fig.  5  Elastic  problem. 

The  solution  to  this  problem,  a.  .(x,t)  and  u.(x,t)  is  known,  and  with 
it,  all  field  quantities  can  be  computed  on  the  boundary.  This  solution 
allows  us  to  define  the  equivalent  elastic  problem  shown  in  Fig.  6. 


u  =  B(x,t) 

n    *  *  ' 


u  = 
n 

e 

u  = 


ue  =0 
s 


u   = 
n 

e 

u  = 


0 
C(x,t) 


an  =  k(E,v)f(x,t) 


Fig.  6  Equivalent  elastic  problem. 

The  transform  of  this  problem,  which  can  now  be  computed  because  the  boundary 
conditions  of  the  moving  crack  are  known  for  all  time,  is  shown  in  Fig.  7. 
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-* 

* 

u       =  B 

n 

* 
s 

ue     =  0 

n 

*         •* 
u      =  A 
s 

ae     =0 

s 

* 
n 

* 
Gij    (X'P) 

* 
u®   (x,p) 

•* 

* 
ct®  =  k(E,v)f* 

u   =  0 


u   =  C 


Fig.  7  Transform  of  equivalent  elastic  problem. 

Now  we  define  an  equivalent  viscoelastic  problem  whose  boundary 
conditions  are  defined  .in  a  manner  to  "correspond"  with  Fig.  7.  This  is 
shown  in  Fig.  8. 


v 


u 


=  0 


v 


ug  -  A(x,t) 


u  =  B(x,t) 

n      '  ' 


a.  .(x,t) 


uV  =  0 
s 


^(Xjt) 


uV  =  0 


u  =  C(x,t) 


V 

ov  - 


o- 


Fig.  8  Equivalent  viscoelastic  problem. 


In  Fig.  8,  K(t)  is  defined  as 


K(t)  =  L-^i 


-ark  (e*,u*)' 


(17) 
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The  transform  of  this  problem  can  also  be  computed  by  virtue  of  the 
definition  of  boundary  values  for  all  time.  The  results  are  shown  in  Fig.  9. 


v* 
u   =  0 


v* 
u   =  A* 


V* 

u       =  B* 
n 

V* 

u       =  0 
s 

V* 

--  0 

V*/                     X 

aij(x'p) 

V*/                     N 

u±  (x,p) 

0" 

s 

V* 

CTs      = 

--  0 

a       = 
n 

=  k(E*,u*)f* 

u. 


V* 

n 


=  0 


v* 
u   =  C* 


an  =0 

Fig.  9  Transform  of  equivalent  viscoelastic  problem. 

Comparing  the  governing  equations  and  boundary  conditions  for  the  problems 
shown  in  Fig.  7  and  Fig.  9?  it  is  seen  that  a  correspondence  exists. 


v*,   v   r  e*l 

ff«(x'p)  =  LaidJs 


E* 


(18) 


u  -*  u1 


This  will  also  be  a  correspondence  principle  for  the  original  problem 
of  Fig.  k   if  it  can  be  shown  that  on  the  extension  in  front  of  the  crack 
the  boundary  condition  of 


o- 


(19) 


produces  u  =0  over  the  same  surface.  It  goes  without  saying  that  the  zero 
shear  stress  condition  must  exist  due  to  the  assumed  symmetry  of  loading  and 
geometry. 
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Using  the  extended  correspondence  principle,  the  displacements  on 
the  boundary  are 


ujT(x,p)  =  [Uf(x'p)TE  -  E*  (20) 


u  -»  u* 

e  e* 

by  symmetry  u  =0,  and  u   =0.   If  there  are  no  material  constants  in 

°     J  J      n  n 

e*  /   n 

u   ,  then  from  (20) , 

n  ' 

u^  =  uS*  (21) 

n     n  v  ' 


and  on  the  boundary 


v    e  /   » 

u   =  u =  0  (22 

n     n  v   ' 


Thus,  the  problem  shown  in  Fig.  8  is  the  same  as  the  original  one  in  Fig.  k 
for  which  a  solution  is  being  sought,  and  a  correspondence  principle  has 
been  establ  jhed  for  this  case. 

Applying  Equation  (l8)  to  the  elastic  stress  solution  of  the  long 
strip  with  displacement  boundary  conditions ,  we  obtain 

v*  (lh     rvo(t)  fiit9(t)]l*l 


E  ->  E* 


Inverting  the  transform, 


k     t       „-v»f,<[e(T)] 


-J  =  ^  f.W^  £  V'-CITJ   V  ^ 


where 


REL 


If  we  look  at  a  point  on  the  line  of  propagation, 

9(t)  =  n  h(t  -  t±)  (25) 

lU 


where  t,  is  the  time  of  arrival  of  the  crack  at  the  point  of  interest. 
Using  Equation  (25),  it  can  be  shown  that  f   in  Equation  (2k)   reduces  to 

f yy[e(t)]  =  h(t)   -  h(t   -  tx)  (26) 

For  a  step  displacement  input,  v  (t)  =  V  h(t),  the  expression  for  opening 
mode  stress  simplifies  to, 

v<-*>  ■  *„  4fjm.  (*  -  o  fe  [  h^(;h.(Te(;g] *   (-) 

Carrying  out  the  differentiation,  and  simplifying 

-I         ♦!       v      A  rEBEL(t)     .  ,t.  W*  '  T)  °(t)    ■     I 
a      (x,o,t)    =  V    thsr  ■-  +      1  dT  (2o) 

yyV  °  V3bTrL^X7TT5    Jo     [2(x  -  c(t)]3/2  J 

With  this  expression  for  hoop  stress,  the  stress   intensity  factor  can 
be  obtained  via  the  definition  of  Equation   (5). 

■*<*) .  x  ^(t)  j*^m  vo  ^jjf&r +  £  i^imv^]  (29) 


It  goes  without  saying  that  this  equation  is  difficult  to  use  in  data 
reduction  for  crack  propagation  tests.  This  has  led  to  a  search  for  an 
approximate  solution  possessing  a  simpler,  more  usable  form. 
APPROXIMATE  VISCOELASTIC  STRESS  DISTRIBUTION 

The  expressions  for  the  elastic  stress  distribution  in  the  neighborhood 
of  a  crack  tip  can  be  written  as 

Ke   (E,a,c) 
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where  superscript  e  denotes  elastic  and  subscript  s  denotes  stationary  crack. 
The  functional  dependence  of  K  upon  E,a,  and  c,  indicates  that  the  stress 
intensity  factor  is  in  general  a  function  of  material,  loading  and  crack 
length.   Because  the  crack  is  stationary,  the  viscoelastic  correspondence 
principle  is  applicable. 

K  f  (E*,a*,c) 
°±]s  =  : t±A*)  (3D 

Inverting  the  transform, 

oj.      =  -=L    (E(t),o(t),c)   f      (9)  (32) 

US       ^2r  ij 

where   as  before,   to  fit  the  definition, 

KiVs  =  L_1  [Kif<E*     *•*>]  (33) 

For  slowly  moving  cracks,  the  inertia  terms  in  the  equilibrium  equations 
can  be  neglected,  and  the  field  equations  for  the  moving  crack  become 
identical  to  those  for  the  stationary  crack.   If  it  were  not  for  the  changing 
boundary  conditions  at  the  crack  tip,  the  1  i  solutions  would  be  identical; 
however  the  S.I.F.  is  altered  as  the  crack  propagates  because  the  geometry  of 
the  body  is  changing  and 


■£  K  <3"o 


In  the  case  of  the  long  strip  laboratory  sample,  the  specimen  approxi- 
mates on  infinite  strip  containing  a  semi-infinite  crack.  As  the  crack 
propagates,  the  geometry  remains  unchanged  to  the  degree  that  the  specimen 
models  the  infinite  sheet;  in  fact,  any  geometry  for  which  the  crack  propa- 
gation makes  minor  perturbations  on  the  nature  of  the  body,  such  as  infinite 

domains ,  will  be  well  represented  by 
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kf(t)  =  [kv]        ,., 

lv  '    L  1SJC  -  c(t) 


(35) 


To  illustrate  this,  consider  an  infinite  sheet  with  a  central  crack  of 
length  2c. 


Kis  ■  °^ 


By  the  correspondence  principle  and  the  definition  of  Equation  (33)? 


(36) 


Kis  =  L_1  [CT*^]  =  o{t)JT 


(37) 


Using  the  approximation  of  Equation   (35)  5   the  viscoelastic  S.I.F.    for  a 
moving  crack  is   found  to  be 


kJ  =  a(t)7cTt7 


(38) 


In  this  case,  this  is  exactly  the  expression  that  results  from  the  more 
rigorous  extended  correspondence  principle. 
Infinite  Strip  With  Displacement  Boundary  Conditions 

Beginning  with  the  elastic  solution  for  a  stationary  crack  in  a  long 
strip, 


n  e  -  /A.  rr    ° 

ijs   v3hrr  L 


V2(x-c) 


By  the  classical  correspondence  principle, 


F.  .(e) 


(39) 


v*        fk 


E* 


^         3bTT     K/2[x^7J 


v  *  _ 
o 


Fid(e) 


(hO) 


Inverting  the  transform,  still  holding  the  crack  stationary, 


v    /4   ij 


T  F.  .(e)  dv 

*  "    J*  W*-*>TT* 


3hn 


72(x-c)  Jo- 


(M) 
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For  a  step  displacement  loading, 


v  =  V  h(t) 
o    o  x  ' 


Ik     Fii(9) 

Converting  to  a  moving  crack  through  the  approximate  relation  of  Equation  (35) 
9(t)  is  defined  as  before, 


v    Ik    „       EREL(t) 


aii  :=  ^Vo  ,7"      [h(t)  ■  ^-Vl  ^3) 

1J    3bn  o  y2(x-c(t)  1 


The  stress  intensity  factor  is  readily  recognizable  from  this  form  as 


^■4ToW*)  m 


In  order  to  measure  this  quantity  on  each  specimen  directly  which  would 

account  for  any  specimen  to  specimen  variation  in  E   (t),  we  seek  to  express 

KEL 

the  S.I.F.  in  terms  of  stress.  For  a  long  strip  without  a  crack,  the  princi- 
pal elastic  stress  is 


V 
e      E     o 

O"   = 


(1*5) 


y   (l-u£)  b 

Applying  the  correspondence  principle  to  an  incompressible  strip, 

<,**  =^E*v*  (k6) 

y   3b    o 


Upon  inversion 


)   +  dv 
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For  a  step  strain  history,  v  =  V  h(t) 


#*>  "  5  vo  W*>  (W 


From  Equation   (1&)  , 


«J<*W£a'(t)  (W 


y 

By  using  the  stress  expression  for  an  uncracked  strip,  a   (t)  represents 
the  stress  in  the  undisturbed  field  far  out  in  front  of  the  crack.   In  order 
to  reduce  laboratory  data,  it  is  necessary  to  relate  the  total  load  on  the 
specimen,  or  the  displacement  of  the  specimen  grips,  to  the  far  field  stress, 
This  has  been  accomplished  through  a  finite  element  program. 

Because  the  stress,  <j   (t),  is  constantly  relaxing  during  the  test,  the 
stress  intensity  factor  is  likewise  decreasing.  This  results  in  the  crack 
velocity  steadily  diminishing  and  approaching  zero  quite  rapidly.  The 
constantly  varying  crack  velocity  and,  more  particularly,  the  fact  that  it  is 
approaching  zero  makes  this  test  somewhat  unattractive;  however  it  is  not 
impossible  to  use. 
Infinite  Strip  with  Stress  Boundary  Conditions 

The  same  method  when  applied  to  the  infinite  strip  with  stress  boundary 
conditions  leads  to  considerably  different  results.  Beginning  with  the 
elastic  solution,  (ll) 


«J-A  <5°> 


The  viscoelastic  stress  intensity  factor  for  a  stationary  crack  can  be 
obtained  from  Equation   (33) 

KV       =  CT(t)  J^  (51) 

IS  V     TT  N        ' 
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Since  the  stress  intensity  factor  is  independent  of  crack  length,  it  is  the 
same  for  a  slowly  moving  crack  as  for  a  stationary  one. 


(52) 


K^(t)  =  CT(t)  J± 


This  agrees  exactly  with  the  rigorous  solution  obtained  in  Equation  (l6). 

Thus  if  stress  is  held  constant,  the  stress  intensity  factor  remains  constant 

and  so  does  the  velocity.   However,  this  test  is  very  difficult  to  run  in 

actual  practice.  Because  the  specimen  is  not  actually  infinite,  the  net 

cross -section  stress  is  difficult  to  hold  constant  as  the  crack  runs. 

(13) 
Hertzler     has  painstakingly  done  this  in  order  to  obtain  data  for  this 

particular  loading  condition;  however  it  is  much  more  convenient  to  work  with 

constant  load. 

Infinite  Strip  With  Constant  Load  Boundary  Conditions 

Probably  the  most  convenient  approach  to  viscoelastic  fracture  character- 
ization from  the  experimental  point  of  view  is  to  apply  a  constant  load.  A 
strip  specimen  is  bonded  between  rigid  grips  as  in  the  case  of  displacement 
boundary  conditions.   (See  Figure  10.) 

1  F  -  Constant 


Fig   IO.      Constant    Force     Sheet 
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The  S.I.F.  for  this  specimen  in  terms  of  the  far-field  stress  is  given  in 
equation  (^9)« 


Ki(t)  =  a(t)  ys  (53) 

The  relationship  between  a(t)  and  the  applied  force  has  been  obtained  from 
a  finite  element  solution  for  a  1"  x  6"  specimen  which  has  been  used  in 
this  research  program.  For  crack  lengths  up  to  three  inches,  the  difference 
between  the  far-field  stress  in  the  undisturbed  region  and  the  average  stress 
based  on  net  cross -section  was  less  than  Vf0.     Considering  this  error  to  be 
negligible , 


^-MM  ») 


where  F(t)  is  the  applied  force  and  A(t)  is  the  net  cross-sectional  area 
computed  from  the  specimen  length  minus  the  crack  length. 

The  variation  of  A  with  crack  length  will  be  very  slight  when  the  crack 
is  small  compared  to  the  sheet  length,  which  should  be  the  domain  of  most, 
if  not  all,  of  the  test.  By  varying  the  level  of  applied  force,  the  velocity 
of  crack  propagation  in  the  sheet  is  varied.  For  a  given  force,  the  S.I.F. 
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is  virtually  constant  as  mentioned  before;  consequently  the  velocity  of 

V        /   \ 

crack  propagation  will  also  be  constant.  A  plot  of  K  vs.  c(t)  completes 

the  characterization. 

Infinite  Strip --Const ant  Displacement  Rate 

A  good  test  case  to  which  the  fracture  characterization  might  be  applied 

is  the  infinite  strip  subjected  to  a  constant  displacement  rate  loading. 

v 
Having  obtained  the  fracture  characterization  of  K,  vs.  6  via  the  constant 

load  test,  it  is  possible  to  predict  the  crack  behavior  for  the  same  strip 

under  constant  displacement  rate  and  check  the  results  experimentally. 

From  the  elastic  solution, 


Using  the  correspondence  relation  of  Equation  (33) , 


1    v3bn   o 
Inverting  the  transform, 

Ik     ,t 
K 


k7*  =  J^r  E*v  *  (56) 


iw-^W^'-ar*  (57) 


For  constant  displacement  rate,  R, 

^-l/i'w*-*  (58) 

o- 

The  integral  can  be  evaluated  by  means  of  a  Prony  Series  representation  of 
the  relaxation  modulus. 

W*>  ■  E  V""1*  +  Er  (59) 
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Substituting  into  the  above 

.E. 


v  AT       r  rhi/  -Qf.tv-i 


(60) 


From  this  expression  and  the  fracture  characterization,  a  crack  history  can 
be  predicted.  This  solution  is  applicable  for  crack  lengths  that  approximate 
the  infinite  strip  condition.   Once  the  crack  becomes  great  enough,  by 
comparison  to  the  length  of  the  strip,  to  cause  the  far-field  stress  to  deviate 
from  the  net  section  stress  when  it  runs,  the  above  expression  will  begin  to  be 
in  error.  The  amount  of  the  error  can  be  determined  numerically. 

FRACTURE  CHARACTERIZATION 

The  first  two  methods  of  fracture  characterization  previously  discussed 

(13) 
have  been  carried  out  experimentally  by  C .  M.  Hertzler    .  The  others  will 

be  carried  out  in  the  near  future. 
Constant  Displacement  Case 

Tests  were  conducted  on  cracked  biaxial  strip  specimens  with  displacement 
boundary  conditions.  The  strips  were  six  inches  in  length,  one  inch  in  width 
(2b)  and  l/lO  inch  thick.  Redwood  tabs  bonded  to  the  horizontal  edges  allowed 

for  gripping.  An  initial  crack  of  0.75  inch  was  cut,  as  suggested  by 

(9) 
Mueller   ,  and  the  length  of  the  propagating  crack  was  measured  as  a  function 

of  time  with  an  optical  comparator. 

The  testing  was  performed  with  the  apparatus,  which  consisted  of  an 
Instron  Universal  Testing  Machine  equipped  with  special  jaws  to  clamp  the 
specimen  and  apply  displacement  v  while  maintaining  the  upper  and  lower 
boundaries  parallel.  The  displacement  was  applied  at  rates  ranging  from  one 
inch/minute  for  low  temperatures  to  five  inches/minute  at  the  high  tempera- 
tures. A  conditioning  chamber  allowed  the  tests  to  be  conducted  at  the  same 
specific  temperatures  as  in  the  previous  case.  After  waiting  ten  times  the 
loading  time  the  ramp  tests  simulated  step  displacement  loading,  v  h(t),  for 
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which  the  stress  intensity  factor  reduces  to  Equation  (kk) .  As  time 
proceeds  after  loading,  the  stress  intensity  factor  is  monotonically 
decreasing  as  a  result  of  relaxation.  This  produces  a  corresponding  decrease 
in  crack  velocity,  which  was  observed  experimentally. 

The  stress  intensity  factor  for  each  test  was  calculated  by  using 
Equation  (U9)  and  plotted  on  a  graph  of  stress  intensity  factor  versus  crack 
velocity  for  each  specific  testing  temperature  (see  Figs.  11  through  16).  A 
conventional  a^,  shift  was  attempted  with  the  data,  and  a  smooth  master  curve 
was  obtained  graphically  by  shifting  the  curves  for  various  temperatures  along 
the  abscissa,  relative  to  the  curve  at  T  to  produce  a  continuous  function. 
The  resulting  characterization  curve  is  shown  in  Fig.  17  where  T  =  70  F, 
and  a_  =  1.   It  turned  out  that  the  a^  curve  for  the  fracture  characteriza- 
tion corresponded  exactly  with  the  a^  curve  for  modulus. 
Consent  Stress 

The  specimen  configuration  used  to  achieve  the  stress  boundary  conditions 
is  shown  in  Fig.  18.   The  specimen  was  cut  from  samples  designed  for  use  in 
biaxial  tests  with  a  gage  length  of  1.0  inch  and  thickness  of  0.1  inch.  The 
wood-  tabs  were  segmented  with  saw  cuts,  which  stopped  in  the  adhesive  bond 
just  short  of  the  propellant .   (To  penetrate  through  the  bond  into  the 
propellant  would  have  created  a  stress  riser  and  a  point  of  probable  crack 
initiation.)  The  width,  w,  ideally  would  be  very  small,  allowing  shear-free 
edges;  however,  practical  considerations  forced  a  width  of  w  =  0.25  inch. 

Initial  experiments,  with  the  right-hand  edge  of  the  propellant  flush 
with  the  wooden  tab  edge,  showed  that  soon  after  loading  a  fast  running  crack 
would  initiate  at  a  local  defect  along  the  right  edge  of  even  a  carefully  cut 
specimen.  Adding  a  l/8  inch  unloaded  length  reduced  the  stress  along  this 
edge  and  eliminated  the  problem.  To  compensate  for  this  departure  from 
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Fig.    18     Specimen  configuration  for  stress  "boundary  conditions 
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uniform  stress  at  the  right  hand  edge  and  yet  still  meet  the  restriction 

leading  to  Equation  (52),  c(t)    was  limited  to  2.0  inches. 

max 

The  thickness  of  the  saw  cuts  in  the  wooden  tabs  was  large  enough  to 
allow  unrestrained  e  ,  or  d/w  ^  e     .  The  upper  segments  were  individually 
suspended  and  the  lower  segments  were  individually  loaded  by  weights  W.  An 
initial  crack  of  one  inch  was  cut  with  a  razor  blade.  Horizontal  extension 
of  the  crack  and  constant  propagation  velocity  were  taken  as  evidence  of 
uniform  loading.  Each  specimen  was  tested  by  applying  all  loads  at  t  =  0, 
and  the  length  of  the  propagating  crack,  c(t),  was  measured  with  an  optical 
comparator.  As  the  crack  passed  a  segment,  the  weight  on  that  segment  was 
removed  in  proportion  to  the  width  of  the  segment  traversed  by  the  crack. 
The  tests  were  conducted  in  an  environmental  chamber  with  the  relative  humidity 
less  than  30$,  at  six  temperatures  ranging  from  15  F  to  85  F.  The  crack 
velocity  was  slower  in  this  geometry  for  a  given  S.I.F.  than  it  was  under 
constant  displacement  as  is  shown  in  Figs.  11  through  16.  The  shifted  master 
curve  for  constant  stress,  which  is  shown  in  Fig.  17,  was  obtained  by  using 
the  same  a—  shift  factor  as  was  determined  in  the  first  set  of  tests  for 
constant  displacement. 

Differences  in  the  fracture  behavior  of  the  load  histories  is  thought  to 
be  due  to  dewetting.  The  theory  developed  considers  only  a  purely 

viscoelastic  material.  However,  the  solid  propellant  was  a  filled  material, 

(1*0 
which  when  loaded,  exhibited  the  additional  mechanism  of  dewetting.  Wood 

has  shown  that  propellant  dewetting  is  not  a  rate  process,  rather  it  occurs 

upon  initial  application  of  stress.  This  behavior  has  a  significant  bearing 

on  the  experimental  results  obtained  in  this  study  as  will  now  be  shown. 

In  both  cases  tested,  the  measured  load  parameter  was  stress,  thus 

automatically  accounting  for  the  modulus  relaxation  in  the  displacement  case. 

However,  the  two  stress  histories  were  quite  different  as  shown  in  Fig.  19. 
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In  the  stress  case,  stress  was  constant  while  in  the  displacement  case  there 
was  an  initial  extremely  high  stress  which  continually  relaxed. 


■*  t 


Constant  Stress 


Constant  Displacement 


Fig.  19  Stress  histories. 

Thus  a  sample  tested  in  the  stress  case  exhibited  minimum  dewetting 
because  of  its  history  of  constant  stress.   However,  in  the  displacement 
case,  for  a  corresponding  stress  at  t  s  10ta  (where  t&  is  the  time  required 
to  apply  the  load)  much  more  dewetting  had  occurred  due  to  the  initial  high 
stress.   Consequently  the  crack  propagated  through  a  material  which  had 
been  significantly  altered  from  that  in  the  stress  case.  The  propagating 
crack  had  fewer  bonds  to  break  per  unit  area,  compared  with  the  stress  case, 
and  accordingly  its  velocity  was  as  much  as  two  orders  of  magnitude  higher. 
This  produced  significantly  different  characterization  curves  for  the  two 
cases  as  shown  in  Fig.  17. 
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Since  dewetting  produces  microscopic  voids  in  the  material,  stiffness 
is  reduced  as  dewetting  increases.  Using  this  fact,  the  increased  dewetting 
of  the  displacement  case  was  verified  in  preliminary  tests  at  15°F.   In 
these  tests,  for  identical  stress  and  time,  the  modulus  in  the  displacement 
case  was  one-half  that  of  the  stress  case. 

In  an  attempt  to  reduce  the  initial  high  stress  peak  in  the  displacement 
case,  thus  reducing  the  dewetting,  additional  preliminary  tests  at  15°F  were 
conducted.  The  maximum  stress  was  limited  to  that  which  produced  stress 
intensity  factors  comparable  with  those  of  the  stress  case  for  the  same 
temperature.  Equivalent  moduli  showed  that  the  amount  of  dewetting  was  in 
fact  equivalent.  However,  as  soon  as  the  displacement  rate  was  stopped,  the 
stress  relaxed  due  to  E_.TTIT(t).  Thus  even  if  the  crack  propagated  at  the  same 
velocity  as  in  the  stress  case  (approximately  0.1  inch  per  hour)  relaxation 
of  the  driving  stress  decelerated  the  crack  to  zero  velocity  over  distances 
too  small  to  be  accurately  measured.  Thus,  no  additional  characterization 
data  were  obtainable  with  this  loading  method. 

The  accuracy  of  the  measured  velocity  in  the  stress  case  was  considered 
superior  to  that  of  the  displacement  case,  especially  for  low  velocities. 
For  in  the  displacement  case,  c(t)  had  to  be  measured  over  relatively  short 
periods  of  time  to  obtain  an  accurate  crack  propagation  function.   (Velocity 
was  obtained  by  graphically  measuring  the  slope  of  this  function.)  As  c(t) 
became  small,  Ac(t)  for  the  time  interval  mentioned  above  also  became  small. 
Thus,  because  the  experimental  crack  propagation  was  irregular  due  to 
inhomogeneities  in  the  propellant,  the  actual  velocity  was  obscured  as  those 
irregularities  approached  the  magnitude  of  Ac(t).  However,  in  the  stress 
case,  c(t)  was  a  constant  and  Ac(t)  could  be  measured  over  large  time  intervals. 

The  scatter  in  the  characterization  curves  is  due  to  the  inaccuracies 
mentioned  above  for  the  displacement  case  and  also  to  the  inhomogeneous  nature 
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of  a  highly  filled  viscoelastic  material  such  as  propellant.  By  compari- 
son, similar  characterization  of  solithane,  an  unfilled  and  relatively 

(9) 

more  homogeneous  material,  by  Francis  et.  al.    showed  significantly  less 

* 
scatter. 

Referring  to  the  characterization  curves,  in  both  cases  there  appears 

to  be  a  K,  critical  below  which  a  crack  will  not  propagate.  Again  for 

both  cases,  especially  apparent  in  the  stress  case,  there  appears  to  be  a 

maximum  velocity  of  propagation;  however  these  maximum  velocities  are 

somewhat  misleading.   Preliminary  tests  showed  that  at  very  high  loads  the 

failure  was  general  and  catastropic,  initiated  at  voids  and  other  stress 

risers  throughout  the  specimen  as  well  as  at  the  crack  tip.   Clean  crack 

propagation  was  possible  at  high  loads  for  low  temperature,  but  high  loads 

at  room  temperature  precipitated  fracture  throughout  the  specimen 

resembling  disintegration. 


Solithane  at  room  temperature  has  a  relaxation  time  of  less  than  one 
minute.  Thus,  though  tested  for  the  displacement  case,  essentially  constant 
velocity  crack  propagation  in  these  tests  allowed  the  additional  advantage 
of  increased  accuracy  of  the  measured  crack  propagation  function  for  reasons 
mentioned  above. 
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ROCKET  MOTOR  ANALYSIS 

The  purpose  of  all  of  the  preceding  material  relative  to  making  a 
fracture  characterization  is  to  obtain  the  base  information  required  to 
analyze  a  rocket  motor.  Rocket  motor  geometries  are  generally  complex 
and  finite  element  methods  are  required  to  generate  solutions  and  stress 
intensity  factors  even  for  elastic  fracture  predictions.  A  method  will 
now  be  presented  which  uses  the  elastic  methods  to  generate  the  desired 
viscoelastic  fracture  results. 

Viscoelastic  fracture  analyses  of  rocket  motor  geometries  compounds 
the  two  complexities  of  time  dependent  stress  analysis  and  fracture 
mechanics.  The  expense  of  the  undertaking,  as  well  as  the  difficulty  of 
doing  it  rigorously,  suggests  that  approximate  methods  of  solution  be 
pursued  and  developed. 

Appendix  I  has  been  written  to  give  motivation  to  an  approximate 
method  of  stress  analysis  commonly  referred  to  as  quasi-viscoelastic.  In 
this  development,  a  circular  cylinder  is  used  as  the  motor  geometry,  and 
it  is  assumed  that  the  method  holds  for  all  geometries.  With  a  means 
available  for  obtaining  the  stress  field,  another  approximate  analysis 
procedure  is  now  presented  to  determine  fracture  behavior. 

For  most,  if  not  all,  two-dimensional  geometries  with  through-the- 
thickness  cracks,  the  stress  intensity  factor  has  a  common  form. 


L  =  Xo-  JT  (61) 


where  X  is  a  modifying  factor  containing  information  relative  to  geometry, 
and  in  a  sense,  reflects  the  departure  of  the  geometry  in  question  from 
the  infinite  sheet  with  a  central  crack  where  X  =  1.  X  includes 
influences  such  as  boundaries,  stiffeners,  dimensionality  of  the  stress 
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field,  etc.  a   is  the  applied  uniaxial  stress,  and  when  cast  in  this  form 
it  is  seen  that  a   represents  the  stress  that  would  exist  in  the  body  at 
the  location  of  the  crack  if  there  were  no  crack  present.   In  other  words, 
a   is  the  stress  acting  across  the  line  of  the  crack  in  a  identical 
geometry  but  without  a  crack.  For  the  rocket  motor  shown  in  Fig.  20, 


Fig.  20 

the  crack  extends  along  the  r  direction,  and  the  stress  acting  across  the 
crack  is  a  .  The  stress  intensity  factor  for  the  motor  can  be  written  as 


K 


1  =  \a     J~c 


(62) 


where  the  crack  length  is  c,  and  the  influence  of  the  bore,  case,  trans- 
verse stress  and  material  properties  are  included  in  X  .   For  a  given 
configuration;  the  functional  dependence  of  the  quantities  in  Eq.  (62)  can 
be  represented  in  terms  of  these  parameters  which  will  vary  with  time  in  a 
viscoelastic  fracture  analysis. 


Kle  =  X(E,c)  ae(E,p)  JT 


(63) 
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The  subscript,  e,  designates  elastic  stress  intensity  factor,  while  E,  c 
and  p  represent  the  modulus,  crack  length  and  pressure  loading  respectively. 
These  are  the  quantities  that  will  vary  with  time  in  a  viscoelastic  fracture 
analysis. 

It  is  proposed  that  \(E,c)  be  found  by  numerical  methods  using  finite 

(2) 
element  codes.  Such  a  study  was  conducted  by  Francis   ,  et.  al.  for 

elastic  materials  and  they  experimentally  confirmed  predictions  made  in 

this  manner.  For  the  geometry  specified,  the  elastic  stress  intensity 

factor  was  computed  from  energy  differences  for  runs  with  incrementally 

different  crack  lengths. 


K' 


EAU 


le   rrt  A  c 


(6*0 


(The  procedure  is  defined  in  a  note  by  Deverall  and  Lindsey    .)  Results 

are  shown  in  Fig.  21  for  a  high  modulus  material.  Similar  curves  could  be 

constructed  for  the  range  of  modulus  values  exhibited  by  propellant,  and  X 

follows  by  dividing  K.  by  a     reevaluated  at  r  =  a  +  c. 

-Le     y 


0.25      0.50      0.75      1.00 
CRACK  LENGTH,  IN. 


Fig.  21 
Stress  Intensity  Factor  vs.  Crack  Length 
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Two  examples  for  low  modulus  materials  were  available  from  the 
literature  in  reference  2  for  pressure  loads.  They  are  reproduced  in 
Figs.  22a  and  22b ,  and  from  them  \   was  computed  to  show  the  dependence  of 
this  factor  on  modulus .   Of  course  this  range  of  modulus  is  much  greater 
than  that  found  in  rocket  motors.  However  once  these  results  are  obtained 
\    is  easily  computed. 

A  plot  of  the  results  in  Fig.  23  shows  that  \    is  not  strongly  sensi- 
tive to  changes  in  modulus.  Two  more  runs  in  the  range  of  modulus 
indicated  by  the  dotted  lines  would  provide  the  information  needed  to 
construct  the  viscoelastic  solution  that  will  now  be  developed  in  detail. 
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Consider  first  a  stationary  crack  in  this  geometry.  We  have  shown 
previously  that  the  viscoelastic  stress  intensity  factor  for  this  case 
may  be  obtained  from  the  correspondence  principle  (Eqn.  33). 

-1 


KnV  (t)  =  L"1  [X*(E*,c)  a>(E*,p*)  /c"]  (65) 


v, 


a  At)    can  be  found  from  the  elastic  solution  by  the  method  outlined  in 
H 

Appendix  I. 

The  viscoelastic  correction  factor,  \    ,  can  be  found  from  Fig.  23 
by  using  an  equivalent  elastic  modulus,  which  can  be  found  for  the  viscoelastic 
problem.   To  illustrate,  the  equivalent  modulus  for  plane  strain  is  defined 
as  follows : 


v=llu[(l.u)CTV_  ^  (66a) 


eq 


In  other  words,  the  instantaneous  viscoelastic  stress  and  strain  can  be 
used  to  find  the  value  of  Young's  modulus  which  would  produce  those 
stresses  and  strains  if  the  body  were  elastic.  Of  course  the  equivalent 
modulus  varies  with  time. 

The  stresses  and  strains  generally  vary  also  from  point  to  point 
within  the  motor;  however  for  a  cylinder  the  sum  of  normal  stresses  and 
strains  is  a  constant  throughout  the  motor.  We  will  therefore  choose  to 
define  the  equivalent  modulus  in  terms  of  this  sum.  Adding  Eqns .  (66a) 
and  ( 66b ) , 


v         v        (l+u)(l-2u)   _   v         v,  /£„n 

er  +  ee  = e [CTr  +  CTe]  (67) 

eq 


k2 


The  viscoelastic  strains  can  be  found  in  terms  of  viscoelastic  stresses 
through  the  viscoelastic  constitutive  law. 


S  +  <  =  (l+u)(l-2u)  J  Dcrp  (t-r)  ±   («£  +  <£)*       (68) 

0-  ' 


Eliminating   (e^  +  6g)   fromEqns.    (67)  and   (68),  the  equivalent  modulus 
becomes, 


v+  v 

Eeq(t) -=-2b (69) 

J  Bcrp  (t-r)A  (a%a> 

0- 


With  this  parameter  determined,  Fig.  23  can  be  used  to  find  \    (t)  =  \[E      (t),cl 

v        eq   '  J 

An  inversion  of  Eqn.  (65)  for  a  fixed  crack  length  yields 

t 
K^(t)  =  JT  J  \v(t-r)  a^(r)dT  (70) 

0- 


where 


\v(t-x)  =  0(E,c)]  E  =  Eeq(t-r) 


Referring  back  to  Eqn.  (69), 


aV(t-T)  +  a!"(t-T) 

v*-t)  ■  — — r       fl    ™ 


With  the  geometry  and  material  correction  factor  included  in  the 
analysis  so  that  a  continuous  updating  of  these  properties  can  be  effected, 


^3 


the  correlation  between  static  and  dynamic  viscoelastic  stress  intensity- 
factor  would  follow  from  Eqn.  (38). 


K^(t)=7cTtTJ  0(E,c)]c=c(t)    aJ(r)dT 

(t-x) 


(72) 


E=E 


eq 


Schematically  the  solution  may  be  viewed  to  proceed  as  shown  in  Fig. 
2h.     The  static  viscoelastic  solution  is  given  for  a  variety  of  crack 
lengths.  The  viscoelastic  stress  intensity  factor  for  a  moving  crack  is 
assembled  from  the  incremental  pieces  of  the  static  solution. 
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Fig.  2k 
Schmetic  of  Solution 
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Procedure 

To  illustrate  how  the  solution  should  proceed,  a  step  by  step  set  of 
instructions  will  be  given. 

1.  An  appropriate  time  increment  is  selected  and  a  and  a„  are  computed 

r     9 

for  p  =  p(t  )  and  E  =  E_.„T(t  )  from  Eqn.  (1-19)  for  any  given 
pressure-time  trace.   (r  is  taken  to  be  a  +  c_.)  A  time  incre- 
ment is  selected  by  guess  using  estimates  of  load  and  resulting 
crack  velocities  as  guidelines.  The  time  increment  should  allow 
for  moderate  extension  of  the  crack  but  not  be  so  large  that  the 
accuracy  of  the  solution  is  compromised. 

2.  Form  a  +  aQ  and  compute  the  equivalent  modulus  from  Eqn.  (69) 

r    0 

for  t  =  t.  and  t  =  t  .   In  summation  form  it  becomes, 

CTr<V  +CTe(tn> 
E      (t    )   = 


eq     n'  n 

V/ ,      \    .       V/  ,      \  s  V/ 


2     D        (t    -t.)[(o-V(t.)+a:(t.))  -a(t.    _)  -a:(t.-l)] 
_     crpv   n     i/LV   rv    iy        91  rv    l-l'        9      1       J 


i=o  crp 


3.  Knowing  E   's  for  the  first  time  increment,  \  's  can  be  found  for 
&  eq  v 

the  corresponding  times  from  Fig.  23. 
k.     Knowing  \      from  step  3  and  aa   from  step  1,  compute  Kn  (t  )  from 

V  0  -LS   -L 

Eqn.  (70)  for  the  first  time  increment. 

1=0 

5.  Using  the  fracture  characterization  curve  for  the  grain  propellant 
(for  example,  see  Fig.  17)  find  the  crack  velocity  associated  with 
the  first  time  increment,  c  .  Knowing  the  velocity  and  the  length 
of  the  time  increment,  a  new  crack  length,  c^,  can  be  computed  to 
be  used  for  the  next  iteration. 


J+5 


Vi 


V/ 


7. 
8. 

9. 
10, 

li, 


Progress  to  t  =  t_  and  c  =  c  and  find  a   (t  )  and  a   (t  ) ,  using 

P  =  P('tp)j  E  =  Erft^?^  and"  r  =  a  +  c"l  for  the  second-  "time 
increment  and  values  as  before  for  t  . 

Compute  E   for  t  =  t_.  t  =  t.,  -  t  =  t_. 

eq.         2'      1'      0 

Find  corresponding  \    's. 

Compute  K  (t  )  and  find  the  new  crack  velocity  and  crack  length, 

The  process  is  repeated  until  the  entire  crack  history  is 

recorded  as  shown  in  Fig.  2k. 

From  this  information  a  crack  history  plot  can  be  constructed  as 

shown  in  Fig.  25. 
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Fig.  25    Viscoelastic  Crack  History 
This  approximate  solution  is  obtainable  with  realistic  computer  times 
and  associated  expenditures.  At  the  same  time  it  should  yield  good  results 
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SUMMARY  AND  CONCLUSIONS 

The  "basic  definitions  for  viscoelastic  stress  intensity  factors  are 
set  forth  with  associated  computations  for  determining  them  in  laboratory- 
specimens.  Laboratory  tests  have  been  conducted  to  characterize  propellants 
fracturewise,  and  the  most  suitable  test  found  was  the  long  strip  subjected 
to  a  constant  force  loading.  The  influence  of  dewetting  on  the  fracture 
characterization  appears  to  be  significant,  but  the  matter  has  not  yet  been 
fully  resolved. 

An  approximate  method  for  performing  viscoelastic  fracture  analyses 
has  been  presented  which  incorporates  the  fore -mentioned  fracture  character- 
ization.  It  requires  minimum  computer  usage  and  should  yield  good  results. 
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APPENDIX  I 

An  approximate  method  of  viscoelastic  analysis  can  be  developed 
from  linear  superposition  of  elastic  solutions.   Such  quasiviscoelastic 
solutions  can  give  good  solutions  if  properly  used.   The  guidelines  for 
use  are  formulated  in  this  appendix  by  developing  the  quasiviscoelastic 
equations. 

The  elastic  normal  stress  expressions  for  an  encased  circular 
cylinder  are  of  the  form 

CTr,9e  =A1P'  +B1P  (I-1} 

where  p'  is  the  interface  pressure,  and  p  is  the  bore  pressure,  A.  and 
B,  are  functions  of  motor  geometry  and  the  radial  coordinate;  however 
they  are  not  time  dependent. 


For  ar  , 


Al         ^2  2   U    2  J 

b  -  a       r 


a2        b2 

Bl  =  ~T       2  [1  "  ~2  ] 
b  -  a      r 


For  cj  , 


i.2         2 

Al  =  ~2 2  t1  +  "5  1 

b  -  a      r 


a2       b2 

Bl  =  ~2 2  [1  +  ~2  ] 

b  -  a      r 


From  the  elastic  analysis, 


p  "  A2  +  B2E 


(1-2) 


where  Ap  and  B?  are  composed  of  geometrical  constants,  and  E  is  Young's 
Modulus  of  the  propellant.   The  normal  stress  expression  can  then  be 
written  as 


ar,9e  =TTm     +  B1P 


(1-3) 


1+8 


By  the  correspondence  principle,  the  transformed  stress  for  a  unit 
step  bore  pressure  is 

r,8l    A  +  BE*   s  U  H; 

where  the  asterisk  indicates  the  Laplace  Transform  and  p*  =  l/s.  The 
subscript  I  on  stress  indicates  the  indicial  stress  response;  i.e.  the 
stress  resulting  from  a  unit  step  input  loading. 

For  convenience  in  manipulation  the  first  term  in  equation  (l-2) 
will  be  designated  o"j  . 

ci*  =  TT^bW  (I"5) 

Clearing  Equation  (l-5)  of  fractions  and  inverting  to  the  time 
plane, 

t  do- 

AOi(t)  +  BCTl(°)  EREL(t)  +  B  L    EREL(t  "  T)  "d^.dT  "  1  (I"6) 

If  the  integral  is  broken  up  into  time  increments  for  evaluation  by  sum- 
mation, the  expression  (l-5)  becomes 

Aaz(t)  +  B»j(o)  EREL(t)  ♦  B  jx  Erel  (t  -  t±)  Aan  =  i  (i-7) 
When  t  =  t  =0,  Equation  (l-7)  yields 


o 


CTi<°)  -  a  ;  b\<o)  ^ 


REL 
This  is  the  same  as  Equation  (1-3)  with  E  replaced  by  EREL('t)' 


h9 


Proceeding  on  in  time  to  t  =  t  ,  Equation  (1-7)  now  becomes, 

A°I(V  +  Bai(0)  W*l>  +  B  EREL(0)iai(V  "  1  <**> 

where 

Aa(tl)  =  ai(t1)  -  ai(o) 

If  BO\.(t.)  E   (t.)  is  added  and  subtracted  to  Equation  (l-9),  it  can 
be  written  in  the  form 

[A  +  B  EREL(tl)]  W  "  *  WVl^l  "  X        (I"10) 
where 

A  WV  "  WV  "  EeeL(°) 
If  A  E_,_T  and  A  o      are  chosen  to  be  small  compared  to  E__T  and  o  in 
the  numerical  computing  process, 

A  EEEL(tl)  =  A  +  B  Ej^V  (1-11) 

Once  again  this  of  the  form  of  Equation  (1-3)  with  E  replaced  by  Er)„T(t1). 

KEL   1 

This  procedure  can  be  repeated  for  t  =  t_;  however  more  manipulating  is 
required.   The  result  is 

[A  +  B  WV1   W   -   B[i  WV   '  W   +  "  W*2)4  °I<V 

+  A  E^tg)   A0l(t2)]   =  1  (1-12) 

By  the   same   criterion  as  before, 


W   =     A  +  B  I        It,)  (I"13) 


RELV    2 
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It  is  assumed  that  this  process  can  be  repeated  for  an  arbitrary 
number  of  time  increments.  As  a  result,  a  general  relationship  can  be 
deduced 

From  equation  (I-U)  the  viscoelastic  stress  components  become  upon 
inversion 


ar,ei(V  s  A  ;  B  4L(tt)  +  Bi         (I"15) 


Proceeding  back  one  step  further,  the  viscoelastic  stress  field  resulting 
from  a  unit  step  in  bore  pressure  is  found  from  equation  (1-3). 


Since  o"  Q  is  the  indicial  response,  the  viscoelastic  stress  for 
an  arbitrary  pressure  history  can  be  obtained  via  the  linear  superposition 
integral. 

v  ,    „    r>t   v 


.;J(t)=j^arV)ei(t-T)MliaT  (l_17) 

In  terms  suited  to  numerical  computation, 


n   v 


ar?e(t)  =  p(o)arj9l(t)  +  £  -^(t-VAp.      (1-18) 
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Relating  this  back  to  the  elastic  stress  expression  through  equation 
(1-16) , 

The  viscoelastic  stress  is  obtained  by  subdividing  the  time  into  increments, 
determining  the  associated  increments  in  applied  pressure,  Apj_  ,  and  then 
appropriately  summing. 


52 


REFERENCES 

1.  E.  C.  Francis,  G.  H.  Lindsey  and  R.  R.  Parmerter,  "Pressurized  Crack 
Behavior  in  Two-Dimensional  Rocket  Motor  Geometries,"  J.  Spacecraft 
and  Rockets,  Vol.  9,  No.  6,  Jun.  1972,  p.  1*15. 

2.  E.  C.  Francis,  et.  al. ,  "Application  of  Fracture  Mechanics  to 
Predicting  Failures  in  Solid  Propellants ,"  AFRPL-TR-70-105,  Air  Force 
Rocket  Propulsion  Laboratory,  Edwards,  California,  Sep.  1970 

3.  W.  G.  Knauss,  "Stable  and  Unstable  Crack  Growth  in  Viscoelastic  Media," 
Trans.  Society  of  Rheology,  Vol.  13,  No.  3,  1969,  p.  291. 

h.     M.  L.  Williams,  "On  the  Stress  Distribution  at  the  Base  of  a 
Stationary  Crack,"  J.  Applied  Mechanics,  Mar.  1957,  p.  109. 

5.  H.  M.  Westergaard,  "Bearing  Pressures  and  Cracks,"  J.  Applied  Mechanics, 
Vol.  6,  Trans.  A.S.M.E.,  Vol.  6l,  1939,  p.  A-U9-A-53- 

6.  G.  R.  Irwin,  "Analysis  of  Stresses  and  Strains  Near  the  End  of  a  Crack 
Traversing  a  Plate,"  J.  Applied  Mechanics,  Sep.  1957,  p.  36l. 

7.  G.  A.  C.  Graham,  "Two  Extending  Crack  Problems  in  Viscoelasticity 
Theory,"  Quarterly  of  Applied  Mathematics,  Vol.  27,  No.  k,   Jun.  1970, 

p.  ^97. 

8.  W.  G.  Knauss,  "Delayed  Failure—The  Griffith  Problem  for  Linearly 
Viscoelastic  Materials,"  Int.  J.  of  Fracture  Mechanics,  Vol.  6,  No.  1, 
Mar.  1970,  p.  7. 

9.  H.  K.  Meuller,  "Stress  Intensity  Factor  and  Crack  Opening  for  a  Linearly 
Viscoelastic  Strip  with  a  Slowly  Propagating  Central  Crack,"  Int.  J.  of 
Fracture  Mechanics,  Vol.  7,  No.  2,  Jun.  1971,  p.  129- 

10.  H.  K.  Mueller  and  W.  G.  Knauss,  "Crack  Propagation  in  a  Linearly 
Viscoelastic  Strip,"  J.  Applied  Mechanics,  Jun.  1971,  p.  *+83. 

11.  G.  A.  C.  Graham,  "The  Correspondence  Principle  of  Linear  Viscoelasticity 
Theory  for  Mixed  Boundary  Velue  Problems  Involving  Time -Dependent 

53 


Boundary  Regions,"  Quarterly  of  Applied  Mathematics,  Vol.  26,  1968, 
p.  167. 

12.  W.  G.  Knauss,  "Stresses  in  an  Infinite  Strip  Containing  a  Semi- 
Infinite  Crack,"  J.  Applied  Mechanics. 

13.  CM.  Hertzler,  "Viscoelastic  Fracture  Characterization  of  a  Solid 
Propellant , "  Naval  Postgraduate  School  Thesis,  Jun.  1972. 

Ik.      G.  H.  Lindsey  and  J.  E.  Wood,  "An  Isotropic  Theory  for  Dewettable 

Solids,"  Naval  Postgraduate  School  Report  NPS-57Li71011A,  Jan.  1971. 

15.   L.  I.  Deverall  and  G.  H.  Lindsey,  "A  Comparison  of  Numerical  Methods 
for  Determining  Stress  Intensity  Factors,"  J.  Basic  Engineering, 
Jun.  1972,  p.  508. 


5k 


DISTRIBUTION  LIST 


No.  of  Copies 


1.  Arnold  Adicoff  10 

Michaels on  Laboratory 
Naval  Weapons  Center 
China  Lake,  California  93555 

2   Defense  Documentation  Center  20 

Cameron  Station 
Alexandria,  Virginia  2231^ 

3 •   Library  2 

Naval  Postgraduate  School 
Monterey,  California  939^0 

h.     Department  of  Aeronautics  1 

Naval  Postgraduate  School 
Monterey,  California  939^0 

5.  Superintendent  1 
Naval  Postgraduate  School 

Monterey,  California  939^0 

6.  G.  H.  Lindsey  5 
Associate  Professor 

Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  939^0 

7.  LT  R.  G.  Kish  1 
Aero  Engineering  Programs 

Naval  Postgraduate  School 
Monterey,  California  939^0 

8.  LT  R.  E.  Kapernick  1 
Aero  Engineering  Programs 

Naval  Postgraduate  School 
Monterey,  California  939^0 

9.  Eugene  Francis  1 
United  Technology  Center 

1050  East  Arquez  Avenue 
Sunnyvale,  California  9U086 

10.  Harold  Zitzer  1 

United  Technology  Center 
1050  East  Arquez  Avenue 
Sunnyvale,  California  9U086 


55 


11.  Dr.  Lamar  Deverall 
United  Technology  Center 
1050  East  Arquez  Avenue 
Sunnyvale,  California  9U086 

12.  Dr.  S.  A.  Murch 
United  Technology  Center 
1050  East  Arquez  Avenue 
Sunnyvale,  California  9U086 

13.  Dr.  John  M.  Wozencraft 
Dean  of  Research,  Code  023 
Naval  Postgraduate  School 
Monterey,  California  93940 


56 


UNCLASSIFIED 


St'.    UMt\      C'l.l  ^>»1  lli'.ll  li<Il 


DOCUMENT  CONTROL  DATA  -R&D 


rrttv  >  ;.is  wf/<  tit  inn  ol  title,   hotly  of  tths  trar  t  and  indexing  ■inriofvif/on  tnttst  hi    entered  when  tin-  itvt-rnll  fpnri  r  .  <  f.i      tli-  ./, 


(1,1  ...4N  a  1 1  Nii    a  c  i  i  v  i  r  y  (Corporate  author ) 

Naval  Postgraduate  School 
Monterey,  California     939^0 


da.   HLf'ONT    SECURITY    CLAtSI 

Unclassified 


2b.     Gf-tOUP 


Rl    I'OHT     TITLE 


Viscoelastic  Analysis  Method  for  Rocket  Motors  Containing  Cracks 


i     DESCRIPTIVE   NOTES  (Type  ol  report  and,inclusive  dates) 

Technical  Report         NP3-57LI73031A 


AUTHORIS)  (First  name,  middle  initial,   last  name) 

Gerald  H.  Lindsey 


REPORT  DATE 

March  1973 


7a.     TOTAL    NO.    OF    PAGES 


.58 


7b.    NO      OF    RE  ES 

15 


a.    CONTRACT    OR    GRANT    NO. 


b.    PROJEC  T   NO. 


9a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.    OTHER    REPORT   NO(S)  (Any  other  numbers    that  may  be  assigned 
this  report) 


0.    DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  Unlimited. 


I.    SUPPLEMENTARY    NOTES 


12.    SPONSORING    Ml  LI  TAR  Y    ACTIVITY 

Naval  Weapons  Center 

China  Lake,  California  93555 


3.  ABSTRACT 


The  basic  definitions  for  viscoelastic  stress  intensity  factors  are  set  forth  with 
associated  computations  for  determining  them  in  laboratory  specimens.  Laboratory  tests 
lave  been  conducted  to  characterize  propellants  fracturewise ,  and  the  most  suitable 
-est  found  was  the  long  strip  subjected  to  a  constant  force  loading.  The  influence 
)f  dewetting  on  the  fracture  characterization  appears  to  be  significant,  but  the  matter 
las  not  yet  been  fully  resolved. 

An  approximate  method  for  performing  viscoelastic  fracture  analyses  has  been 
rresented  which  incorporates  the  fore -mentioned  fracture  characterization.  It  requires 
linimum  computer  usage  and  should  yield  good  results. 


)D  .FN°.?..1473    (PAGE  " 

/N  01 01 -807-681 1 


UNCLASSIFIED 


57 


Security  Classification 


A-31408 


UNCLASSIFIED 


Si.un:\   Classification 


key    wo  ro  s 


Viscoelastic 

Fracture 
|  Solid  propellant 
;  Fracture  characterization 

Crack  propagation 


LINK      A 


-INK        C 


DD  .'.I™  .1473  'back 


UNCLASSIFIED 


-  /-'      U  1  <i  1   -  10  7  -  6  H  ?  1 


58 


Security  Classification 


DUDLEY  KNOX  LIBRARY 


3  2768  00391410  2 


I     . 


